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Abstract 

We present a computational LCHO-CI approach allowing for the simulation of exchange inter- 
actions in gated lateral quantum dot networks. The approach is based on single-particle states 
calculated using a linear combination of harmonic orbitals (LCHO) of each of the dots, and a con- 
figuration interaction (CI) approach to the interacting electron problem. The LCHO-CI method 
is applied to a network of three quantum dots with one electron spin per dot, and a Heisenberg 
spin Hamiltonian is derived. The manipulation of spin states of a three-spin molecule by applying 
bias to one of the dots is demonstrated and related to the bias dependence of effective exchange 
interaction parameters. 

PACS numbers: 73.21.La,73.23.Hk 
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I. INTRODUCTION 



The electron spin is a quantum two level system, a natural candidate for a qubit.-'^*^ 
For this reason there is currently significant interest in coupling of individual electron spins 
localized at different spatial locations to realize quantum gates.-i^i^ The coupling of the 
spins of localized interacting electrons is investigated using the Heisenberg Hamiltonian.- 
In this approach, the Coulomb interactions between electrons are reduced to the exchange 
coupling of their spins. The exchange interaction is parametrized by the exchange constants 
J dependent upon the specific implementation of the system. 

In lateral gated quantum dot devices the spatial localization of individual electrons can 
be achieved by electrostatic coupling to their charge. To date, controlled confinement of 
one, two, and three spatially separated electrons has been demonstrated using the single,-"^ 
double and triple^ii^ quantum dots, respectively. In these devices the electronic 
orbital degrees of freedom are manipulated directly by tuning the gate voltages, which, 
only through topology and statistics, translates into the control over the total spin of the 
system.— 'i^ii^'i^ii^'i^ The voltage control of exchange coupling of two electrons in a single 
dot^S and coherent control of spin states of a two-electron double-dot molecule^^ was re- 
cently demonstrated. The Heisenberg Hamiltonian appropriate for this two-spin system is 
parametrized by a single exchange constant J, whose magnitude defines the energy gap be- 
tween the spin singlet and triplet eigenstates. By measuring this gap one can establish the 
dependence of J on the gate voltages. In a more complex network of quantum dots with 
one spin per dot the number of exchange couphngs needed is equal to the number N{N — l)/2 
of pairs. Even the simplest network of quantum dots, a triple quantum dot with one electron 
per dot, is described by a Heisenberg Hamiltonian with three exchange constants, depend- 
ing nontrivially on the geometry of the system and the gate layout. This dependence was 
studied, e.g., by Scarola and Das Sarma,— who used the Hubbard, variational, and exact 
diagonalization approaches to demonstrate that the three-spin model is valid only for a lim- 
ited range of triple-dot parameters. Mizel and Lidar— ^i^ii^ arrived at similar conclusions 
using the Heitler-London and Hund-Miilliken schemes to calculate the energy levels of three 
coupled dots with one electron per dot. In both cases the many-body effects were responsi- 
ble for the appearance of higher-order terms in the effective spin Hamiltonian. The above 
studies were performed on triple-dot systems consisting of three dots on resonance, and did 
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not account for tuning of individual dots. In Ref. 



13 we used the Hubbard model, and in 



Ref. 



171 - the real-space wave functions coupled with the configuration-interaction technique 



(RSP-CI) to analyze the voltage-tunable three-electron gated lateral triple-dot device, how- 
ever without mapping the resulting electronic properties onto the three-spin model. The 
Hubbard model is simple but requires parametrization of the Hamiltonian and cannot be 
directly hnked to gates and applied voltages. The RSP-CI technique is very accurate but is 
difficult to implement for the future simulation of the time evolution of the quantum system. 

The purpose of this work is to present a different computational approach, where linear 
combination of harmonic orbitals (LCHO) of each of the dots is used to describe tunneling 
of electrons, and the configuration interaction (CI) approach is used for the treatment of 
exchange and correlation in the interacting electron problem. The LCHO-CI method can 
then be directly used for the derivation of the effective Heisenberg Hamiltonian consistent 
with gate voltages, and for the simulation of the quantum operations. 

We illustrate our method by analyzing the tunability of the exchange interaction and 
manipulation of the three spin system with voltage in a triangular quantum dot molecule. 
The triangular quantum dot molecule, realized recently by Gaudreau et al.,^^ is needed both 
for implementation of the quantum teleportation and creation of three-particle maximally 
entangled GHZ state. 



II. LCHO-CI METHOD FOR THE ELECTRONIC STRUCTURE OF QUANTUM 
DOT NETWORKS 



A. Single electron states in a quantum dot network 



We describe here the LCHO method for the calculation of single electron states in a 
network of two-dimensional quantum dots. The method is illustrated on the example of a 
triple quantum dot. 

We consider a triple quantum dot molecule created electrostatically by lateral gates^^ and 
approximate its lateral confinement potential by a sum of three two-dimensional Gaussians, 



with Xi and yi being the coordinates of the center 



of ith Gaussian (i=l,2,3), Vi being its depth, and di being its characteristic width. The 
depth and width of each Gaussian are functions of the gate voltages. The centers (xj, 
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of each dot are arbitrary but in the rest of the paper they are chosen to he in the corners 
of an equilateral triangle. In the following we express all distances in effective Bohr radii 
a*Q = efi^ /rrfe^ and all energies in units of effective Rydberg Ry* = e^/2eal, where e and 
m* are the electron charge and effective mass, respectively, and e is the dielectric constant 
of the material. GaAs parameters m* = 0.067 tuq and e = 12.4 give Ry* = 5.93 meV and 
a*Q = 97.9 A. Figure [1] shows the triple-dot confining potential for Vi = V2 = Vs = 10 Ry* 
and di = d2 = = 2.3 a*^. The dimensionless Hamiltonian for one electron in the potential 
of the triple quantum dot network is written in the following form: 



2 3 



dx'^ dy"^ ^ 6 p 



(x - Xj)^ + {y- yif 
d? 



(1) 



i=l 

Since an exact analytical solution of the eigenvalue problem of the Hamiltonian ([T]) is not 
known, an approximate method is needed to calculate the energies and eigenvectors of one 
electron in the quantum dot network. Here we employ a quantum dot analog of the Linear 
Combination of Atomic Orbitals (LCAO) method. Expanding the Gaussian potential of 



dot i to second order in r, — ^exp — I'-* ^ — \^ + Vi , results in the harmonic 
oscillator (HO) potential —Vi + ^Qf{r — riY, with Qi = 2y/Vi/di. The eigenfunctions of this 
HO potential, (t)]^^{x,y) = fn{x — Xi)ipl^{y — yi), are products of the ID HO eigenfunctions 
Mt) = (fe)' - |)"exp , where I = h = = (|)'. The first three 

states have the exphcit form (poit) = (^) " exp (^—^^ , = (^) "* (j) exp ; 

Hence, the molecular eigenfunctions of Eq. ([1]), are written as linear combinations 
of the HO orbitals (LCHO) centered on each dot: 

3no 
i=\ 

where are the expansion coefficients and no is the number of HO orbitals per dot. To 
simplify notation, the three indices n, and m of the HO eigenfunctions (pn^ni^'V) have 
been replaced by the composite index i. Substituting this expression into the eigenvalue 
problem of Eq. ([T]) and multiplying on the left by {(pj \ gives 

3no 3no 



\H\<p,)a^ = enJ2(<i',m<- 



1=1 i=l 
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Defining as the Hamiltonian matrix in LCHO basis with elements (^qoI-^I^oo) 
J dr(f)Ql{x,y)H(f)lfQ{x,y) and 5*,^ as the overlap matrix with elements (^qoI^oo) 
J dr(f)l*Q{x,y)(f)QQ{x,y) allows us to write the generalized eigenvalue problem as 



(3) 



Now, defining a new vector 



6" = (V^^)cr (4) 

and multiplying Eq. ([3]) by (a/S^)^^ on the left gives the standard eigenvalue problem 

(5) 



In order to calculate {^/S^)~'^ the eigenvalue problem of the overlap matrix S(i,Vs = VsEs 
is solved. Here Vs is the matrix with the eigenvectors and Es is the diagonal matrix with 
the eigenvalues. Once the values of these two matrices are calculated, is obtained 

from {^yS^)~^ = VsE^^^'^Vg . Then the energies e„ ([3]) and eigenvectors ([2]) of the electron 
in the quantum dot network can be calculated using Eqs. ([5]) and (jlj). The accuracy of the 
solution depends on the number Suo of HO orbitals included in the basis. Increasing the 
number of basis states increases the accuracy of results. 

We now analyze the matrix elements of the Hamiltonian in the basis 

{\^oq)^ \^oo)^ I0oo)} showing exphcitly the various contributions to the diagonal onsite en- 
ergies and off-diagonal tunneling matrix elements: 



Ha 



^1 tl2 tl3 
^21 ^2 ^^23 



^31 ^32 



(6) 



If we label each of the three dots with indices i,j,k, we can express the diagonal matrix 
elements as 



^00 



H\ 



'00 



+ a) + {<Plo\SV^\4>lo) + {<Pm\V'\4>lo) + (0ool^l0oo), 



(7) 



where we denote —Vi exp 



V^Q + 6V^, with V^Q being the 2D harmonic 



oscillator potential associated with dot i: V^q = —Vi + ^Qf{r — fi)"^. The first term in 
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©, {—Vi + Qi), is the dominant term. It is the energy of the ground state of the har- 
monic oscillator potential V^q. The second term, (0ool'^^*l</'oo)' gives the correction due to 
non-harmonicity of the confining potential. The third and fourth terms, (0ool^"' l</'oo) 
(^qqI^^I^qq), give the correction to the energy level of an isolated quantum dot i due to the 
presence of the other two quantum dot potentials. 

The off-diagonal matrix elements describe electron tunneling from dot i to dot j. The 
tunneling matrix elements are determined by several contributions : 

= i-V, + m'PMm) + (0'oo|5^i0oo) + (0'ool^'l0oo) + (0'ool^l0oo)- (8) 

The first term is directly proportional to the overlap between HO wave functions centered on 
different dots, (0001*^00)' ^^"^ therefore small for dots which are far apart. The second term, 
(000 1 I ^00)' is correction due to non-harmonicity of the confining potential. The third 
term, (0ool^"' I'/'oo)) i^ ^ "two-centered" integral, as it involves a product of three functions 
centered at two different dots. The fourth term, (0ool^'^l*/'oo)' i^ ^ "three-centered" integral 
being thus the smallest term. This analysis of matrix elements can also be applied to a 
larger HO basis including more shells {p, d, . . .). 



B. Many-electron states in a quantum dot network 

The many-electron Hamiltonian of the quantum dot network written in second quantiza- 
tion is 

^ = Y1 ^i^^i^i + \ l](^^>I^Oc^c+CfcQ, (9) 

where and Cm are, respectively, the creation and annihilation operators of a particle on 
the spin-orbital itinerant (molecular) state m, ipmif) = ^ni^^Xms, where Xms is the electronic 
spinor corresponding to the spin = ±1/2. Indices i,j, k and / run over spin-orbitals, ej is 
the single-particle energy of an electron in the molecular state j and v is the dimensionless 
Coulomb potential, ^ ^ , . 

f ; \r2-ri\ 

The Coulomb matrix elements (CMEs) in the itinerant basis are computed as linear 
combinations of CMEs in the localized HO basis 

Srio S'lio 3no Srto 

{ij\v\kl) = {Xi\Xi){Xj\Xk)^^^^alaiatal{rs\v\tu), (10) 

r=l s=l t=l u=l 



where {rs\v\tu) = J dn J df2(j);{ri)(j)*{r2)j^^^\M^2)(puiri). Using the identity = 
^ / Y exp [iq{fi — r2)], we have 

1 /»oo /'27r 

{rs\v\tu) = - dq ddFXiq.mtuMKtMGStM (H) 

Jo Jo 

with F^^{q,d) = J^^dtipl{t)ip„,{t) exp [±itq cos 1?] and G^^(g,^9) 

(ity^* exp [±ztg sin These four integrals were obtained analytically and 

the double integral (ITT]) over q and -t? was carried out numerically. 

In order to obtain the energy levels and coherent eigenfunctions of electrons in the 
quantum dot network we combine the LCHO method with the configuration interaction (CI) 
approach. In the CI method we build the iV-electron basis out of all possible configurations 
of the N electrons on the molecular single-particle states. The many-electron eigenvalues and 
eigenfunctions are obtained by diagonalizing the Hamiltonian (Q in this A^-electron basis 
of configurations. The number of these configurations riQ (and hence the dimension of H) 
depends on the number of electrons and the number nso of molecular spin-orbital states 
through He = i^nso-'Ny.Ni • ^^^^ as = 3 electrons, LCHO-CI calculations 

using only s-type HO orbital per dot gives 6 spin-orbitals, leading to ric = 20. If for LCHO 
calculation we also include p-type HO orbitals, there are 18 molecular spin-orbitals, giving 
ric = 816, while including also d HO orbitals gives ric = 7140. The number of configurations 
in the 3-electron basis increases very rapidly with increasing number of HO orbitals in the 
one-electron calculations thus increasing the computational requirements. 

One is primarily interested in quantum networks with one electron (one spin) per quantum 
dot. If one retains only one orbital per dot, the number of spin-orbitals is nso = 2A^, and 
number of possible configurations is reduced to nc = -^^tm- -^^^ = 10 spins we already 
have 184756 configurations. 

Further reduction in the number of configurations is possible. Since the Hamiltonian 
is rotationally invariant, the Hamiltonian needs to be diagonalized only in one of the 
subspaces of lowest \Sz\. In what follows we present numerical results for the three-electron 
Sz = —1/2 and 5*^ = —3/2 states obtained with the LCHO-CI method. 
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III. EFFECTIVE HEISENBERG MODEL 



In order to establish a connection with quantum computation, it is convenient to approx- 
imate the Hamihonian of the system of singly-occupied dots in a quantum dot network 
by an effective Heisenberg spin Hamiltonian. For three electrons in three quantum dots this 
Hamiltonian takes the form 



1^^ 1^^ 1^^ 



(12) 



where Jij are the dimensionless exchange coupling constants, and a are the Pauli matrices. 
From this Hamiltonian the quantum gate V SWAP can be obtained by turning on one Jij 
for an appropriate amount of time.- Combining this gate with single qubit gates we can 
obtain the CNOT gate, and CNOT together with one-qubit gates form the universal basis 
for quantum computing.- 

Since commutes with Sz-, matrix is block diagonal in the basis of eigenvectors of Sz- 
Here we will treat the subspace Sz = —1/2. Taking the basis vectors {| lit)? I iti)?! Tii)} 
we obtain 



1 
4 



J 



12 



-J23- 
2^23 
2Jl3 



J 
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2J23 2Ji3 
— J12 — J23 + <^13 2Ji2 
2Jl2 —^12 + J23 



Jl 



13 



(13) 



Let us introduce a Jacobi basis of spin states: \(3a) = -^{\ lit) - I iti)), IPb) = -^(1 lit 
) + I iti) - 2| tii)), and 1/3,) = ^(| iit) + | iti) + I Tii)) • The Heisenberg Hamiltonian 
(IT^ in the Jacobi basis has the following form: 



~^Jav — |(^23 — Jav) "^(^12 — Jis) 



^(Jl2 — J13) 










4 •-'ai; 



(14) 



where Jav = {J12 + J23 + Jis)/^ is the average exchange constant. We see that if all exchange 
constants are equal the two states \Pa) and \Pb) are degenerate eigenstates with energy — | Jqu 
and total spin S = 1/2, while the state |/3c), corresponding to total spin S = 3/2, has energy 
+ ^Jav The Hamiltonian of the two degenerate states \Pa) and \Pb) is analogous to the 
Hamiltonian of a single spin in a magnetic field, where the magnetic field in the z direction 
is proportional to 7(0/23 — Jav) and the magnetic field in the x direction is proportional to 



4 ("^12 — -'^is)- This is why our system of three electrons in three dots can be thought of 
as a single coded qubit, whose logical states can be manipulated by varying the exchange 
constants.— 

IV. RESULTS FOR A TRIPLE QUANTUM DOT NETWORK 

We illustrate our theoretical approach on an example of a triple dot molecule with pa- 
rameters given in Fig. [TJ In section IIV Al we present the energies and eigenf unctions of 
one electron in the triple dot molecule calculated with LCHO method. In section IIVBI the 
energies of three electrons in the triple dot molecule obtained with LCHO-CI method are 
analyzed and compared with results obtained with Hubbard model.-'* In section IIV CI we 
analyze the triple dot molecule with one and three electrons when one of the dots is being 
biased. Finally, in section IIVDI we use the Heisenberg model to study the effects of biasing 
of one dot in a triple dot molecule with three electrons and compare the results with those 
obtained using LCHO-CI method. 

A. One electron in a triple quantum dot 

Figure [2] shows the set of one-electron energy levels in the triple quantum dot potential 
of Fig. [T] obtained by solving Eq. for different LCHO basis sets: s : {\(PIq), \4>oo)^ I0oo)}' 

s+P ■■ {|0oo), I0lo), \<Pli)^ I0oo), I0?o), \<Pli)^ I0oo), I0?o), I0oi)}, s+p + d, and s + p + d + f. The 
number of orbitals per shell in the spectra of a triple quantum dot is three times the number 
of orbitals per shell in the spectra of one quantum dot. Thus, we have three s orbitals, six 
p orbitals, nine d orbitals and so on. 

The spectrum of one electron in the HO potential associated with one of the dots has 
electronic shells with level spacing Qi = l^fVijdi = 2.75 Ry*. For the electron in an isolated 
Gaussian dot and assuming LCHO basis sets composed of sp, spd and spdf HO shells one 
obtains s-p level spacings of 2.44, 2.46 and 2.37 Ry* respectively. When the electron moves 
in the potential of the triple dot, the tunneling hybridizes shells from different dots and 
leads to molecular levels. Comparing the spectra of one electron in the triple-dot molecule 
for increasing number of HO shells in the LCHO basis we find that as the number of shells 
increases the low energy molecular levels decrease their energy and converge to a definite 
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value. The ground state of one electron in the symmetric triple-dot molecule is nondegen- 
erate, and the first excited state is doubly degenerate. Figure [3] shows the eigenfunctions of 
these three levels (s orbitals ^i(^),^2(^) and C,3{r)) calculated using HO basis sets with dif- 
ferent number of orbitals. Results obtained from calculations with only s-type HO orbitals 
in the LCHO basis are qualitatively the same as those obtained when more HO orbitals {p, 
d, f) are included. Based on this conclusion, and anticipating a larger number of electrons 
in a quantum dot network, we will build the configurations for the three electron prob- 
lem including one-electron energies and eigenfunctions obtained from calculations with the 
s-LCHO basis. Here we show typical numerical values for this s-LCHO calculation. The 
off-diagonal overlap matrix elements appearing in Eq. Q are s = (0ool*i^oo) ~ 0.004088. 
The diagonal on-site energies defined in Eq. ([7]) are e'^ = —8.647493 Ry* and have the con- 
tributions {i ^ J ^ k): {-Vi + VLi) = -7.250193 Ry*, (0^ol'^^'l</'oo) = -0.166187 Ry*, 
(0j,o|1/^|0^o) = -0.615556 Ry*, and (0^ol'^''l0oo) = -0.615556 Ry*. On the other hand, 
the tunneling matrix elements hybridizing the atomic orbitals of adjacent dots, defined in 
Eq. (JED, are t = -0.067162 Ry* and are composed of {i ^ j 7^ k): {-Vi + r],)(</'ool</'oo) = 
-0.029641 Ry*, (0^ool'^^i0oo) = -0.014140 Ry*, (0^ool'^-' I0oo) = -0.018489 Ry* and 
(</'ool^'^l</'oo) — —0.004892 Ry*. After solving Eq. ([3]), three energy levels corresponding 
to hybridized s shell and their eigenfunctions are obtained. The s-shell energy gap ob- 
tained between the ground state, and the degenerate excited state, 1^2) and 1^3), is 
Aes = 0.095 Ry*. We plot these energies in Fig. [2] (first spectrum on the left). The eigen- 
functions of these three levels, plotted in the left-hand part of Fig. [3l are: 

lei) = l/v/3(l + 2s)(|0i) + 102) + 103)), 

16) = l/yW^)i\h) - 102)), (15) 
les) = l/v/6(l-s)(|0i) + |02)-2|03)), 

where s is the off-diagonal overlap matrix element. This number appears because the basis 
eigenfunctions are not orthogonal. 

B. Three electrons in a triple quantum dot network 

Once the one-electron problem has been solved, we proceed to solve the three-electron 
Hamiltonian ([9]) using the configuration interaction method (CI). The first step is to build 
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the basis of three-electron configurations in which H is to be diagonahzed. This is done 
for the subspaces Sz = —1/2 and Sz = —3/2. There are nine three-electron configurations 
with Sz = —1/2 obtained by distributing the three electrons among the six molecular spin- 
orbitals i (T)))"^ = 1?2,3} (for illustration see Fig. 111(a)). These configurations can 

be grouped into six doubly-occupied configurations: 1^4) = Cj'j^c^j^c^lO), \B) = c^^c'^^c'^^\0) , 
|C) = c^]4|C+ |0), \D) = C+C+C+ |0), \E) = C+C+C+ |0), and \F) = c+cjc+ |0), and three 
singly-occupied configurations: \G) = c^|Cj]"|C^-|.|0), \H) = c^C3J^c^^|0), and |/) = c^c^^cj|||0). 
For the spin-polarized system with Sz = —3/2 there is only one possible configuration, 

\K) = 4,4,4,10). 

The construction of the Hamiltonian matrix in these bases of configurations re- 
quires the knowledge of Coulomb matrix elements. The Coulomb matrix elements in the 
itinerant basis are computed from Coulomb matrix elements in the localized basis as de- 
fined in Eq. (ITO!) . using the coefficients of the one-electron molecular eigenfunctions writ- 
ten in HO basis. The localized CMEs, Eq. ( fTTl) . for the triple quantum dot potential of 
Fig. [1] consist of the onsite repulsion U = {rr\v\rr) = 2.939179 Ry*, interdot repul- 
sion V = {rs\v\sr) = 0.512857 Ry*, and a number of small electron-electron scattering 
terms: {rr\v\rs) = 0.004653 Ry*, {rs\v\st) = 0.003446 Ry*, {rs\v\rs) = 0.000049 Ry*, 
and {rr\v\st) = 0.000019 Ry*. Upon rotation into the itinerant basis we find that the 
largest Coulomb matrix elements include direct repulsion on the lowest kinetic energy level 
(ll|i;|ll) = 1.315811 Ry* and on the degenerate excited state (22|i;|22) = (33|v|33) = 
1.730884 Ry*. The direct and exchange interaction terms between the ground and excited 
states are (12|t;|21) = 1.320137 Ry* and (12|i;|12) = 0.806983 Ry*, respectively, while the 
scattering term (ll|f |22) = 0.806983 Ry*. Similar matrix elements are obtained for the 
degenerate shell of excited states: (23|t;|32) = 0.918383 Ry*, (23|i;|23) = 0.406251 Ry*, and 
(22|u|33) = 0.406251 Ry*. These Coulomb matrix elements are used in the construction 
of the three-electron Hamiltonian ([9]). After diagonalizing this Hamiltonian in the basis 
of configurations with Sz = -1/2: = {\A),\B),\C),\D),\E),\F),\G),\H),\I)}, we 

obtain a spectrum of nine levels, shown in Fig. [5] (left). The nine eigenfunctions are 
linear combinations of the basis configurations = Yli=A^i\'^)- spectrum we can 

clearly distinguish two groups of levels separated by a large gap. The group in the upper 
part of the spectrum is composed of six energy levels and the group with lowest energy has 
three levels. We focus on this low-energy group. It is composed of a doubly degenerate 
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ground state and a non-degenerate first excited state, as shown in Fig. H] (b), separated 
by the energy gap AE = 0.0027 Ry*. The doubly degenerate ground state has total spin 
S = 1/2 and its two eigenfunctions \A') and \B') have their biggest contribution from the 
three electron configurations \A) and \B), respectively. However, these states are highly 
correlated, with large contributions from other configurations. The eigenfunction \A'), with 
biggest contribution from the state 1^4), has also contributions from configurations \F),\H) 
and |/): \A') = 0.6014|A) + 0.5533|F) - 0A069\H) - 0.4082|/), while the eigenfunction \B') 
with biggest contribution from state \B), has also contributions from states \C), \D) and 
\E): \B') = 0.6014|5)-0.4075|C)-0.5533|D) + 0.4075|E), as shown in Fig. 11(c). The first 
excited state of this Sz = —1/2 three-electron spectrum, state \GHI'), has total S = 3/2. 
Its eigenfunction \GHI') has approximately equal contribution from the three singly occu- 
pied molecular configurations: \G), \H) and |/). Taking into account all four Sz subspaces, 
the ground state is four times degenerate (two levels with Sz = —1/2 and two levels with 
Sz = 1/2) and the first excited state is also four times degenerate (one level belonging to 
each Sz subspace). 

In order to compare the microscopic model with the Hubbard model^^ we propose a mod- 
ified LCHO-CI method. First, we carry out the s-LCHO one-electron calculations neglect- 
ing the overlap matrix. The resulting one-electron spectrum also shows a non-degenerate 
ground state and a doubly degenerate first excited state but with a bigger energy gap 
Aef"^ = 3t = 0.201 Ry*. Next, we apply the CI method to the Sz = -1/2 subspace 
of Hamiltonian ([9]) but neglect all localized CME terms different than U and V. As in 
the itinerant electron basis, the three-electron spectrum obtained in this calculation shows 
two groups of levels, shown in the right-hand part of Fig. [51 The lower group conserves 
the structure of levels (a non-degenerate ground state and doubly-degenerate first excited 
state), although the energy gap is bigger AE^"^^ = 0.0111 Ry*. On the other hand the upper 
group of levels not only shows bigger energy gaps, but the fifth and sixth excited states, 
nondegenerate in the full LCHO-CI calculation, become degenerate. 

C. Triple dot under bias 

In this section we study the evolution of the energy spectrum of a three-electron molecule 
as a function of Vi, i.e., the depth of dot 1. The applied bias is kept smaller than the s-p 
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energy gap of one electron in one quantum dot in order to prevent population of the biased 
dot with two electrons. This is done so that we can attempt to map the three electron 
spectra onto the three spin spectra obtained in the Heisenberg model. We carry out our 
calculations for the case of s-HO orbitals as LCHO basis. Figure M (a) shows the resulting 
one-electron spectra as a function of bias of dot 1. The bias is measured in the units of the 
one-electron energy gap Aes = 0.095 Ry*. When the dot is biased, the degenerate excited 
levels split due to the breaking of symmetry. This energy splitting increases with increasing 
bias . It is also observed that all energies decrease, but each of the three energies at a 
different rate. This can be understood by analyzing the evolution of orbitals associated with 
these energies. When the three dots are on resonance, the two excited levels are degenerate, 
and any linear combination of the two eigenfunctions 1^2) and 1^3), defined in Eq. f ll5l) . is 
also an eigenfunction of the system corresponding to the same energy. This leads to equal 
probability of finding of the electron in each of the dots. However, if one of the dots is 
different, the excited level splits into two, and the wave functions 1^2) and l^s) reflect this 
symmetry breaking. The first excited state 1^2) has a contribution from the HO orbital 
corresponding to the first dot that decreases with increasing bias Vi (at the same time the 
contribution of this orbital increases in the ground state). On the other hand, the second 
excited state l^s) does not have any contribution from this orbital. 

Figure int^b) shows the spectra of three electrons with Sz = —1/2 in the triple dot molecule 
for different values of Vi. We are interested in the three lowest levels as these are the ones 
relevant for quantum information processing and thus the ones that will be mapped onto 
Heisenberg spectra. Figure [71(a) shows the energies of these levels measured from the ground 
state energy. The degenerate ground state splits with bias. When bias is increased, the 
energy gap between the ground state and the first excited state increases, while the energy 
gap between the first and second excited states decreases. 

D. Heisenberg model for three electrons 

We now attempt to map the three-electron spectrum for the case in which dot 1 is biased 
onto the spectrum of the three-spin Heisenberg Hamiltonian (fT4l) in the Jacobi basis. The 
mapping relies on the assumption that biasing dot 1 can be modeled by bias-dependent but 
equal parameters J12 = J13 = J, while J23 is different, but it may also be a function of 
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bias. With Ju equal J13, the two Jacobi states \Pa) and \Pb) are eigenstates, with energies 
-^Jav - 1{J23 ~ Jav) and -f + 1(^23 " Jav), respectively, and \(3c) is an eigenstate with 
energy +jJav The three eigenstates can be written in a way which emphasizes the role of 
the electron in the dot under bias (dot 1): \(3a) = \ l)\S) and \Pb) = 1/V3| l)\To) - T 
)|T_). Here the singlet and triplet states involving the second and third dot were written as 
\S) = 1/V2{\ it) - I TD), \To) = 1/V2{\ it) + I Ti)) and |T_) = | ii). 

Figure [3(b) shows a plot of the energy spectrum of the Heisenberg Hamiltonian as a 
function of J/J23 from J = J23 to J = O.93J23 with J > (antiferromagnetic ground state). 
By changing J from J23 to J = we drive the system from three equal dots to dot 1 totally 
decoupled. We compare this spectrum to that obtained in the full electronic calculation, 
shown in Fig. [7](a). As we can see, the two spectra behave in a similar manner as a function 
of bias of one of the dots ranging from three equal dots (zero bias) to a bias of ~ 1.5 times 
the one-electron energy gap. We can now propose a procedure of finding the parameters 
Jij by fitting the energy gaps in the Heisenberg spectrum, Ae^^ = J23 — J and Ael^ = f"^) 
to those in the LCHO-CI calculation. The value of J23 can be extracted at zero bias from 
the gap between the doubly-degenerate ground and the excited states. Assuming J23 to be 
independent of the bias of dot 1, two values of J - one from each gap of the spectrum - can 
be obtained and averaged for each step of bias. The numerical results show a decreasing 
value of the coupling constant J as the bias is being increased. This is the expected behavior 
because as we bias dot 1 the exchange of the electron in dot 1 with electrons in dots 2 and 
3 decreases. 



V. CONCLUSIONS 



We presented a computational LCHO-CI approach allowing for the simulation of ex- 
change interactions in gated lateral quantum dot networks. The method was illustrated 
by analyzing the electronic properties of a lateral triple quantum dot network with one 
electron per dot. The LCHO-CI calculations show a low-energy spectrum composed of an 
antiferromagnetic {S = 1/2) ground state separated by a small gap from the spin polarized 
{S = 3/2) excited state, and separated by a large gap from the remaining excited levels in- 
volving double occupancy of quantum dots. We have shown that the behavior of these eight 
low-energy levels with bias of one quantum dot can be effectively reproduced by a Heisen- 
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berg spin model for a certain range of bias applied to quantum dots. For each value of bias, 
exchange coupling constants can be calculated from energy gaps in LCHO-CI model. We 
have thus established a connection between physical "external" parameters such as voltages, 
and exchange interaction among spins in the Heisenberg model. 
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FIG. 1: 3D plot of the triple quantum dot potential with equal dots characterized by depths 
Vi = 10 Ry* and widths di = 2.3 a^. The dots are centered at (0,0), (4,0) and (2,3.4641) forming 
an equilateral triangle with side lengths 4 a*^. 
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FIG. 2: Energy levels of one electron in the triple quantum dot potential of Fig. [T] calculated with 
the LCHO method as a function of the number of HO basis states . 
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S - LCHO. - LSHS Wd - LCHO SPBF - LCfciS 

FIG. 3: Comparison of the three lowest energy eigenfunctions IC2) and l^s)) of one electron 

in a triple quantum dot, calculated with the LCHO method using different numbers of HO wave 
functions as basis. Eigenfunctions at the bottom correspond to the ground state, and the two 
at the top, 1^2) and l^s), correspond to the degenerate excited states. 
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S =1/2 3-ELECTRON CONFIGURATIONS: 
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FIG. 4: (a) Examples of three-electron configurations that form the basis for CI calculations, (b) 
The lowest group of levels of the three electron spectrum, (c) Contribution of different configura- 
tions to one of the two states of the ground level. 



20 



-21,8^ 
-21,9- 
-22,0 - 
-22,1 - 

-22,2 - 
-22,3 - 

^ -24.37^ 
111 -24,38- 
-24,39 - 
-24,40 - 
-24,41 - 
-24,42 - 



« 



LCHO-CI 



2x 



2x 



2x 



S=1/2 
S=1/2 
S=1/2 



S=3/2 
S=1/2 



2x 
2x 
2x 



S=1, only V and U 

localized CME 
different than zero 



2x 



FIG. 5: Left: spectrum of three electrons with Sz = —1/2 in the triple quantum dot calculated with 
LCHO-Cl method. Right: spectrum of three electrons with Sz = —1/2 in the triple dot calculated 
with LCHO-CI method but neglecting the overlap matrix and all CMEs except for U = {rr\v\rr) 
and V = {rs\v\sr). 
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FIG. 6: (a) Spectra of one electron in the triple dot as a function of bias of dot 1. (b) Spectra of 
three electrons with Sz = —1/2 in the triple dot a function of bias of dot 1, calculated with the 
LCHO-CI method considering s-HO levels in the LCHO calculation. 
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FIG. 7: (a) Three lowest levels of spectra of three electrons in the triple dot (Fig. [6|), with the 
ground state energy as a reference level, as a function of bias of dot 1. (b) Spectra of three electrons 
Sz = —1/2 given by Heisenberg model with J12 = J13 = J plotted as a function of J/J23 with J 
varying from J23 to O.93J23. 
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